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Abstract. The diagonal spin-spin correlations of the square lattice Ising model, originally expressed as 
Tocplitz determinants, are given by two distinct Frodholm determinants - one with an integral operator 
having an AppcU function kernel and another with a summation operator having a Gauss hyporgcometric 
function kernel. Either determinant allows for a Neumann expansion possessing a natural A-parameter 
generalisation and we prove that both expansions are in fact equal, implying a continuous and a discrete 
representation of the form factors. Our proof employs an extension of the classic study by Geronimo and Case 
1161 . applying scattering theory to orthogonal polynomial systems on the unit circle, to the bi-orthogonal 
situation. 



1. Form Factor expansions of the Diagonal Correlations 

The two-dimensional Ising model and the classical random matrix ensembles are examples of model 
systems with wide ranging applicability. Taken literally, the two-dimensional Ising model specifies a classical 
magnetic system from the microscopic interactions of two-state spin variables. These variables are confined 
to a particular two-dimensional lattice, and the interactions are confined to neighbouring lattice sites. In 
particular, in the case of the square lattice, each site carries a spin variable (Jij £ {—1, 1} with coupling 
between nearest neighbours in the horizontal and vertical directions. The joint probability density function 
for a particular configuration {o-ij}-N<i.j<N of the states on a (2A^ + 1) x {2N + 1) lattice centred at the 
origin is given by 

^ I N-l N N N-l \ 

(1.1) P2N+l{{(T^,J}) = exp /Ci ^ ^ cr,;jCr.,+ij -Hifz X! XI '^''^^'^'J + l ' 



2N+1 



i=~N j = -N i=-N j=-N 



where Z2N+1 is the normalisation such that summing (jl.ll) over all allowed {cij} gives unity. The parameters 
Ki and K2 are dimensionless couplings in the horizontal and vertical directions respectively. As a model of 
a magnet it would be natural to generalise (jl.ip to include a coupling to an external magnetic field which 
adds 

N 

h (Jij 

i,j=-N 

to the exponent in (jl.ip to give the joint probability density P2N+i{Wi-j})- the absence of such a term 
one is said to be considering the zero field Ising model. 

It is a celebrated result due to Peierls that the zero field two-dimensional Ising model exhibits a phase 
transition from a high temperature disordered phase characterised by zero spontaneous magnetisation to a 
low temperature ordered phase with non-zero magnetisation. Setting ^ = sinh 2Ki sinh 2K2 the Kramers- 
Wannier duality argument gives that the phase transition occurs at ^ = 1. The zero field magnetisation is 
specified by 

(1-2) (cro.o) lim lim V cro.o^2W+i({c^«j})- 

Moreover there is a closed form expression for ((To.o) in the low-temperature phase 

(1-3) M = (1 - 

which was announced by Onsager in 1948 [57] and proven by Yang in 1952 [H] . 
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Our interest is in the zero- field spin-spin correlation function between spins cro.o at the centre of the lattice, 
and the spin (T,;,j at site In the infinite lattice limit this is defined as 



(1-4) (cro.oo-jj) = Ihn cro.o<JijP2N+i{{cri.j})- 

N^oc ^ — ^ 

Onsager and Kaufman knew of a Toeplitz determinant form for the case i = j = n (spin-spin correlation 
along the diagonal) but never published their result. However a draft paper has recently come to light which 
obtains this result and, according to Baxter [3], was almost certainly written by the pair. This reads [50] 

where 

1/2 

(1.6) a(0) 



1 - t-'e- 



1 - kr^e'^ 

and n e N. This result (jl.Sp in fact includes (|1.3p as a special case. To sec this one notes that for n large 

upon the assumption that to leading order spins far apart are independent of each other. On the other hand 
for Toeplitz determinants with symbol a(0) such that loga(0) is Laurent expandable in e'^ 

(1.8) loga(0)= ^v^'^'^ 

P— — 00 

and with the further requirement that X]^-oo blkpP < '^'^ Szcgo theorem [3S] gives 

(1.9) det[ai_j]ij=i,...,„ cxp ncg -t- pCpC_p + o(l) . 



For a{&) given by f|1.6p and ^ > 1 (low temperature phase) 

(1.10) co = 0, cp=-^r'^', peZ, p^O. 
The Szego theorem then gives 

(1.11) (ao,o'T„,„) ~ (1 - e')^/^ 

which when substituted in p.7p implies (|1.3p . 

Intimately related to the asymptotic expansion (|1.9p is the so-called form factor expansion 

(1-12) (ao,o^„,„) = (1 - r')^/^ (^1 + 

and its A generalisation 



(1.13) C_(n,n;A)==(l-r')i/M l + ^A^ 



(2p) 
n.n 



Expansions of the type (|1.12p were initiated by Wu [15] in the study of the large n form of p.4p with 
* = Oj J = so the spins are along the same row. Each term in p.l2p is to be a higher order correction in the 
large n limit; and to leading order in is proportional to ^~2p("+p). The coefficients /n^n^ furthermore 
have the interpretation of 2p quasi-particlc contributions to the two-point correlation function [43] 

An explicit 2p-dimensional multiple integral formula for j'i^'^ follows as the special case of the explicit 
form of the form factor expansion of {<7Q_Qan.n) given in [43] . However this is complicated, and in fact starting 
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directly with (jl.Sp the simpler expression 



(1-14) fiV = J^ I dx,--- I dx,. 



Jo 



2p p 



{1 - tX2j-l){xJ_^ - 1) 



P P 

X J]^ ]^(1 - tX2k~lX2jy^ Y\. (^2j-l - 2:2fc-l)^(2^2i - X2k)'^ , 
j=l k=l i<j<k<p 

where t = < 1 and n > 0, has been derived originally by [3], [TO] and independently by [5B], [5]. In the 
high temperature regime the correlations are given by 

oo 

(1-15) Koa„,„) = (l-0'/'E/n'^'^ 

p=0 

and 

Mp+l/2)+p{p+l) rl rl 



2p+l p+1 p 



fc=l j = l ^2j-l ^.=1 
p+1 p 

X ]^(1 - ia;2j-iX2fc)~^ (0:2^-1 - a;2fe-i)^ (x2j-a;2fc)^, 

j=l k=l i<j<k<p+l i<j<k<p 

and in this case t = ^ < 1 and n > 0. 

It is our aim to develop the analytic properties of 12111. 16|) . building on known results from [M], [7|, 
[29j . In particular a theory will be developed which leads to two distinct Fredholm determinant formulas for 
the A generalised form factor expansion (|1.13p . Wc begin in Section [5] by evaluating the first two terms in 
the power scries expansion of (|1.14p and writing (|1.13p as a Fredholm determinant of a particular integral 
operator supported on (0,1). In Section |3] we express (|1.12p as a discrete Fredholm determinant, making 
use of the Borodin-Okounkov identity (16j . [5] from the theory of determinantal point processes fl3j and 
we conjecture that with the kernel multiplied by (|1.13p results. The rest of our study is devoted to 
proving this conjecture and thus the equality of the two Fredholm determinants, although our working is 
more general. Thus our first step in this undertaking is to introduce the bi-orthogonal polynomial system 
attached to a general weight a{C,) in Section 2] Next we indicate in Section [5] how the classic results of 
Geronimo and Case [TO] can be generalised to non-hermitian Toeplitz determinants. Finally in Section IH] 
we pose a A extension of a key functional equation arising in Section |4] and deduce some consequences of 
this, in particular we describe how this is connected with the pioneering work of Wu |42j . It is through this 
connection that wc can establish the desired equality. 



2. Fredholm Determinant with Appell Kernel 

2.1. t ^ Expansions of form factors. Wc see from (im \rm that f'i,l\ f'i,t'^^ have the small t 
expansions 

00 00 

(2.1) f (2P) = iP("+P) f (2P+1) = ^n(p+l/2)+p{p+l) ^>^fc 

fc=0 fc=0 

where c^'^ = c'j^'^{n,p) and t = respectively. Using knowledge of Selberg integral theory (see e.g. [TO] ) 
the first two coefficients of (|2.ip can be readily computed. 

Proposition 2.1. For t = < 1 we have 

< ^ 1 r(n + p + i)r(p + i) ('-^ r(n + j + ^)r(j + i)r(j + 2 ) 

(p!)27r2p r(n+i)r(i) T{n+p + j + l) 
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and 

(2.3) .<=e<J^|4p(„ + p) + ll, 

while for t = ({^ < \ we have 

.34. > ^ 1 r(n + i)r(i) r(n + j + §)r(j + §)r(j + 2 )^ ' 

^ ' ' ^° p!7r2p+i r(n + p+l) 1 JLi T{n + p + j + 2) 



and 

{n + P+\) 
2(n + 2p + l)2 

Proo/. We read oflt from (|1.14p that 



(2.5) =c^^:^^^— ^^-^[4p(p+l)(n + p+i) + n + 2p+l]. 



2p p 



(2-6) co= 



2 /•! /•! ^i-' 

W2Z2^ / ^^2:1 • • • / dxap n n 



(2^27-1 - 1) 



l<j<k<p 



(p\)^7r^P Jo 



dxi I dx. 



4 ■ 




^1 P ^ \ 

k=l l<j<k<p J 



(1-1 P 
/ dxi / dx3--- dx2p^lY\_x'^'^y^'^{l ~ X2k-iy^^^ Yl {X2j-l - X2k~lf 
•^0 -^0 -^0 /c=l l<J<fe<P 

We recognise both integrals in (|2.6p as special cases of the Selberg integral 

(2.7) Sn{\u\2,X) ■■= dh--- dtNY[tf'{l-ti)^- Y[ \t,-tk\^\ 

•'^ 1=1 l<j<k<N 

which, according to Selberg |34j . has the gamma function evaluation 

, ^ ^ _ r(Ai + 1 + J A)r(A2 + 1 + J A)r(l + (j + 1)A) 

Sn{XuX2,x) [[ r(Ai + A2 + 2 + (iv + j-i)A)r(i + A) • 

The result ([Z!2)) follows. In relation to ((231) we read off from ([TTTi)) that 

where, with 5'p(Ai, A2, A)[/] denoting the Selberg integral ^2.7\i with an additional factor / in the integrand 

p 

A, = Spin - i, i, Spin + \. 1), 

A2 = Spin - i, i, l)Spin + i, -i, 1) i,] , 

A3 = 5p(n-i,i,l)[f]t,]5p(n+i,-i,l)[^t,]. 

But according to a result of Aomoto [I] 

1 c M ^ Ai + (p-l)A + l 

^r^^-^5p(Ai,A2,A)[^^f,J -P;,^ + ^^ + 2(p-l)A + 2' 

which gives (|2.3p . The corresponding results for the high temperature phase are found by identical means. □ 
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2.2. Predholm determinants with continuous kernel. Our next result is to show that J2'^=o f^h^i with 
each term as given by (|1.14p can be expressed as a Fredholni determinant with kernel K{x, y) of the so-called 
integrable type [5D] 

(l){x)-ip{y) - 4i{y)ip{x) 



K{x,y) 



x-y 



Proposition 2.2. The low temperature diagonal spin-spin correlation of the square lattice Ising model 
(o'o,ofn,n); \t\ < ^, n>0 is, to within the prefactor il.ll] }. the Fredholm determinant 



(2.8) det(I-K-„,)) = l + ^/(2rt, 
where K~ is the integral operator with kernel K~{x,y) of the explicit form 

(2.9) K-{x,y) = -^^^±Mltf^+\xyr/'^'/' [(1 - x){l - y){l - tx){l - ty)]-'/' 



1 



x-y 



[xFi{n + i; -i, 1; n + 2; t, tx) ~ yF^in + i; -i, 1; 7i + 2; i, ty)\ 



id Fi{a;b,c;d;x,y) is the first Appell function |31j . ^16.13. Consequently we have 



(2.10) 



(2p) 
n.n 



(_1)P /•! 



Jo 



dui ■ 



Proof. Our proof parallels that taken by [3^ in their treatment of a similar integral for the impenetrable 
Bose gas. For convenience we make the replacements X2j-i Uj,X2j n- Vj in (|1.14l) and separate out the 
multiple w-integrations from those of the w-integrations thus 



(2p) 



(p!)27r2p 



dui 



.r^n[(i- 


-tuj){u. 1 


p p 






-to,)(i;-i 



-1/2 



where Ap(u) is the Vandermonde determinant or product of differences of the The inner w-integral 

can be written as a determinant using the Heine identity (see Eqs. (2.27) and (2.2.10) of [35] and the original 
[2T] 1 to give 



(2.11) 



pi det 



/ dv V 
Jo 



n+j+k 



_,[{l-tv){rr^-l)]y^ 



i<j,k<P 



However this determinant form does not yield the simplest evaluation so we employ the additional identity 
for general distinct indeterminates {yj}^^i 



(2.12) 



det 



dx X- 



w{x) 



1 



A2(y) 



det 



dx 



w{x) 



ix~yj)ix-yk) 



l<_7,fc<p 



i<j,fe<p 

Applying this to (j2.1ip we arrive at our simplest form, except for the evaluation of the matrix elements. 
Using the integral representation of the first Appell function, i.e. see [31], Eq. (16.15.1) 



r(Q!)r(7 — a) 



r(7) 



Fi(a; /?,/?'; 7; y) = / duu'^-\l-uy'-'^-\l-xu)-'^ {l-yu)~^' , Re(a) > 0, Re(7-a) > 0, 



with 



Fi{a;(3,l3';r,x,y)= J2 



(«)„+„ (/3)™(/3')n a:" y'^ 



(7) 



m+n 



7. ' T?, ' 



<\x\,\y\)<^, 



we arrive at the result (12.81) and (12.9 



□ 



There is an analogous result for the high temperature correlations although we have in this case the added 
complication of an odd number of integrations and therefore an unpaired variable. 
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Proposition 2.3. The high temperature diagonal spin-spin correlation of the square lattice Ising model 
(o'o,oO'ra,n); 1^1 < 1, n > is, aside from the prefactor a Fredholm minor and its Neumann expansion 

coefficients are 



(2.13) 



(2p+l) 
n.n 



1 

Jo 







/ dvi- ■ 


■ / dvp dct 


h 


Jo 



Ktivk) 



P > 0, /, 



(1) 

n.n 



1 



-K, 



1 



i<j,fe<p 



where the kernels are 

(2.14) K^x, y) = £i!i±|Mlr(x2/)"/^-^/^ [(1 - y){l tx){l ty)f' 

[xFi{n + i; i, 1; n + l;t, tx) ~ yFi{n + i; i, 1; n + 1; i, ty)] 



x~y 



(2.15) 

and 

(2.16) 



if+(x) 



, r(r7, + |)r(|) /2 p ^„ , 1 1 . „ I 1 . -f^ 
Aq — j 1 2-t'i(n + ^, + i,t). 



Proof. For convenience we make the replacements X2j~i '-^ Uj,X2j i— > Vj in (jl.l6p and separate out tlie 
multiple u-intcgrations from those of the u-integrations making the latter the inner ones 



(2p+l) 



J n,n 



^«(p+l/2)-p(p+l) 

p!(p+ l)!7r2p+i 7o 



dv^--- 1 dv, n ^r'""' [(1 - t^^M' - 1)] Ki^) 



X Cdu,--- f du,+ 'f{ u--' [(1 - i^.,)(«7i - 1)]-'/' Al^,iu)Y[ f[iu, - t-\,^') 

•^0 -^0 3 = 1 j=lk=l 

The inner u-integral can be written as a determinant using the Heine identity to give 



(p+ l)!dct 



i<jM<p+i 



We next again employ the general identity (|2.12p but because there are only p factors in the denominator of 
the integrand instead of p + 1 we introduce an extra factor (m — 2)^ into the numerator, which yields 



(2p+l) _ 



J n,n 



^«(p+l/2) 

p!7r2p+i 



/' dv^--- f dv, n v--\l tv,YI\v-^ ~ ^f'^vf^^ n(z;p+i - V,) 
Jo Jo j-^i j-^i 



X det 



and 



dxx''-\l^tx)-^/^ix-^ 



1) 



-1/2. 



J = l 

I{Vj,Vk) I{Vj,Vp+i) 
I{Vp+i,Vk) I{Vp+i,Vp+i) 

{x-t~\l,r 



l<j.k.<.p 



{x - t-'^vj^){x ~ t-^v^^)' 

Naturally the above integral is independent of Wp+i as one can verify by elementary row and column sub- 
tractions. Also it clearly has the structure of a coefficient in the Neumann expansion of a Fredholm minor 
which is not manifestly apparent in the final result (|2.13p . □ 

Remark 2.1. The Fredholm determinant forms found here are implicit in the work of Jimbo and Miwa |23| 
although many details are lacking in this work for us to make a systematic comparison. Suffice it to say that 
we have evaluated the Neumann expansion coefficients of their resolvent kernels, Eqs. (13) and (14), and 
found that they are essentially the same as the ones that can be deduced from (|2.8p and 



Remark 2.2. The Appell functions that appear in our application are the special cases 

Fi(n+ i;i,l;n + l;t,te), Fi(ri + i; -i, 1; n + 2; i, te), 

of Fi{a; P, /3'; 7; x, y) by virtue of relations amongst the parameters /? + 7 = a + 1, /?' = 1, however we are 
not aware of any evidence in the literature |llj . [37j . |38j that they reduce to a single term of univariate 
hypergeometric functions or their products. 
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Remark 2.3. Both expansions, (|1.12p with (|1.14p and (|1.15p with (|1.16p . are vahd for \t\ < 1 and are clearly 

inapplicable at the critical point t ~ 1 as the multiple-integral for each individual form factor fn. diverges 
as i — )■ 1~. The Appcll function appearing in (|2.9|) is well-defined for t = l{'-f — a — f3 = 2) however the 
integrands in (|2.10p are too singular at one endpoint for the integrals to remain finite. In contrast the Appell 
and hypergeometric functions appearing in (|2. 14112. f51) diverge logarithmically (7 — a — /3 = 0) ast— > 1. 

3. Borodin-Okounkov Identity 

For Toeplitz determinants with symbol a{d) such that the Szego formula (jl.9p holds with cq = 0, there is 
a general transformation identity to a Fredholm determinant derived by Borodin and Okounkov [B] in the 
context of determinantal point processes. Later it was realised that the identity had been derived previously 
within Toeplitz determinant theory by Geronimo and Case [16]. In [5] some special cases of the general 
identity were worked out. One of these, appropriately further specialised, is the Toeplitz determinant with 
symbol (|1.6p . Thus from [5] it follows in the low temperature regime n > 

(3.1) (^To,oCT„,„) = (1 - t)'/^ det(I - K„,„+i,...), 
where IK„_„_|_i^... is the integral operator with kernel 

(3.2) K{i,j) 



(-l)m(^),+it(-+^'/^+^ 1 
(l-t) t-j 



^ + i^'^i(-i 5; * + 1; + 2; ^) - -^2Fi(-i J + 1; ^)2^i(i f ; * + 2; ^) 

supported on the successive integers n, n -|- 1, n -|- 2, . . .. Here 2F1 is the standard Gauss hypergeometric 
function and K{i,i) is the limit j — > i of (|3.2p . The Kummcr relation 

2Fi{a,b;c;t) = {I - t)-''2Fi{a,c ~ b;c; j^), 

allows p.2p to be equivalently written 



(3.3) K{l,j)= ^ 2Wl(2)j + l ^(»+j)/2+l 



1 



^^2i^i(i« + |;« + l;02^^i(ij + ^;j + 2;i)--^2Fi(i,.? + |;j + l;02Fi(i,i + i;i + 2;i) 



IJ ■ 

for \t\ < 1. We note that this kernel diverges as t — > 1 because of the fact that one of the Gauss hypergeometric 
functions diverges with the parameter combination c — a — b = —1. 

All Fredholm determinants enjoy natural A generalisations. Thus according to the general Fredholm 
theory |39] we can expand 

00 

(3.4) det(I- A2K„,„+i,...) = 1 + ^(-A2)p ^ dct[K{nj,nk)]i<j^k<p- 

p=l rii>n2>->np>n 

This and (|3.ip compared to (|1.12p and (|1.13p suggest an equality between the coefficients of X^p in the 
expansions (|1.13p and p.4p . 

Proposition 3.1. Let fn,n be specified by {l-H^ and let K{i, j) be given by HS.,'^) . We have 

(3.5) fi%^^i-ir E dct[if(n„n,)]i<,-,<p. 

rii>n2>--->np>n 

We will develop the proof of this statement, which is the main result of our paper, in a series of steps 
starting in the next section. 

A simple piece of supporting evidence in favour of p.Sp is that for p = 1 the leading i — > behaviour of 
the right-hand side agrees with that implied by (|2.ip and (|2.2p . This is seen by noting from (|3.3p that 

^' ^^0 ((z+l)!)2 

and noting furthermore that 

K{ni,ni) ^^Kln^n). 



However, it is not straightforward to extend this analysis to show that p.Sp is consistent with (|2.ip and (|2.2p 
in the t — >■ hmit for p > 1. 

It is known from [21], [H] that 



(3.6) anit) = 



t{t - 1)^ log(ao,oa„,„) - it, T< T, 



t(t- l)^log(ao^oa„,„)-i, T>T, 
satisfies a differential equation which is particular case of the cr-form for Painlcvc VI 

/ , , (i^a \ ^ , ^ da \ ^ da f , da ,\ f da 

subject to the boundary condition 

C_(n, n- A) = (1 - t)V4 + (i|M!^r+i(i + o(t)), 

as t — > 0. In [7], [55] evidence is given that the A- generalisation (|1.13p satisfies the same differential equation, 
with the A-dependence entering only in the boundary condition. 

On the other hand the identity p.ip tells us that logdet(I— K„^„+i^...) satisfies a simple variant of p.7p . 
Proposition [SH] together with the validity of the conjecture from [29] then gives that logdet(I— A^K„^„+i ...) 
satisfies the same differential equation up to the boundary condition. Specifically, this would imply that 
with 

(3.8) a{t) = t(t - 1)4 logdet(I - A2K„,„+i....), 

dt 

the differential equation p.7p is satisfied and furthermore that this equation together with the boundary 
condition 

a{t)^^^-\\n + l) coUit"+\ 

where cq is specified by (|2.2p . completely determines det(I — A^K„,„+i_...). It remains an open problem to 
verify this characterisation. 



4. Bl-ORTHOGONAL POLYNOMIALS ON THE UNIT CIRCLE 

Our first step in the derivation of p.Sp is to build up a theory of bi-orthogonal polynomials on the unit 
circle and its relationship to general, non-hermitian (i.e. complex weight) Tocplitz determinants. The study 
by Gcronimo and Case [16] analysed systems of orthogonal polynomials on the unit circle, or equivalently 
the case of licrmitian Toeplitz determinants and is therefore inadequate for our purposes. However we will 
see that much of their analysis generalises in an obvious manner provided one makes suitable distinctions 
between variables, most notably splitting the orthogonal system into a bi-orthogonal one, and remains valid 
by replacing certain positivity requirements by non-vanishing ones. In his 1966 prescient study of Toeplitz 
determinants applied to the spin correlations of the Ising model Wu [42] uncovered some of the structures 
that exist in systems of bi-orthogonal polynomials without explicitly appreciating or exploiting that fact. 
Much of the theoretical development outlined in this section of our present study draws upon the foundations 
given in fl5| . 

We consider a complex function for our weight w{z), analytic in the cut complex z-plane and which 
possesses a Fourier expansion 

oo 

W{z) = ^ WkZ^, Wk = 
k— — oo 

where z e I? C C and T denotes the unit circle \C,\ = 1 with C = e'^, 9 G (— tt, tt]. Hereafter we will assume 
that z^'w{z),z^w'{z) G L^(T) for all j € Z. We will also assume that the trigonometric sum converges in 
an annulus D = {z E C : Ai < \z\ < A2} and T C D. The doubly infinite sequence {'Wk}^_ao the 
trigonometric moments of the distribution w{e^^)dO /2tt and define the trigonometric moment problem (see 
[35j for an up-to-date account of the herniitian case = W-p). In [16] their variable Z = z = e^^ and their 
measure dp{9) is related to our weight by 

dz 

dp{9) = w{z)-, 
iz 
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and consequently their moments c„ = Wn- In the work of Wu [42j his weight C{z) = w(z) and so his moments 
Cn = Wn- In the Borodin-Okounkov work we identify their weight 4'{z) = w{—z) and V{z) — logw{—z). 
Utihsing the trigonometric moments we define the Toeplitz determinants n > 1 



(4.1) IniCwiO] := det 



(4.2) 



where e will take the integer values 0,±1. In |16| their Toeplitz determinant D„ = /„+i. 

We define a system of bi-orthogonal polynomials {'^n{z),(pn{z)}'^^Q with respect to the weight w{z) on 
the unit circle by the orthogonality relation 

(4-3) / TT-jWiOVniiO'^niO = 5m.n- 



27r< 

This system is taken to be orthonormal and the coefficients in a monomial basis are defined by 



i=0 



(p„(z) = K„Z" + /„Z" ^ + m„z" ^ + . . . + </3„(0) = ^ C„.,Z' 



n.j '^'^ 7 



where k„ is chosen to be equal to k„ without loss of generality (this has the effect of rendering many results 
formally identical to the pre-existing theory of orthogonal polynomials). In |16j the variable K{n) — k„. 
Notwithstanding the notation c„ j in general is not equal to the complex conjugate of Cn.j and is independent 
of it. We also define the reverse or reciprocal polynomial by 

n 
J=0 

The bi-orthogonal polynomials can equivalently be defined up to an overall factor by the orthogonality with 
respect to the monomials 

(4.4) /^'«(C)^n(C)C^ = 0<J<n-l, 
whereas their reciprocal polynomials can be similarly defined by 

(4.5) fi^MCMiO'U = l<J<n. 

The polynomials defined in |16| are related to ours by 

(j){Z,ji) ^ ipn{z), (t>{Z,n) = ifn{z), (j)*{Z,n) = ip*^{z). 

In [12 his Eqs. (2.7), (2.12) are directly comparable with our orthogonality condition (|4.5p . and consequently 
his polynomial X{z) = k„(/?* (z) after reconciling the normalisations. 

The linear system of equations for the coefficients c„ j-, c„j arising from 

(4.6) Cn,n / :r^w(C)¥'«(C)C = r 

27ri(; II m — n 

/ . /" , ^, _ , I < m < 71 — 1 

(4.7) c„,„ / -^w{orvn{o = r - - 

7t ^T^iQ 1 m — n 



has the solution 



/ Wo 



det 



'-rA 



Wn-l 



Wo I 



det 

/ Wo 



( Wo 
\Wn 



Wo 



det 



W-n+l 
\ W-n 



Wi 
Wo J 



-ri] 



( 



det 



Wo 



\W-r 



Wo, 



and in particular one has the following results. 

Proposition 4.1 f P].|15]). The leading and trailing coefficients of the polynomials ipn{z), (^n(z) are 

1 ln[w(0] 



K„ /„+l[w(C)] ' 



CnO = Vn(0) = (-1)^ 



1 /n[Cu;(C)] 



C„0 = <^ri(0) = (-1)'' 



1 In[C-'w{C)] 



K„ /„+l[w(C)] ' Kn /„+l[w(C)] 

The following existence theorem for the bi-orthogonal polynomial system is due to G. Baxter. 

Proposition 4.2 The bi-orthogonal polynomial system {(/3„(z), ^ni'^)}n'=o ^^^sts if and only if In[w{(^)] ^ 

for all n > 1. 

The system is alternatively defined by the sequence of ratios rn = ipn (0) /k„ , known as reflection coefficients 
because of their role in the scattering theory formulation of orthogonal polynomials on the unit circle, together 
with a companion quantity r„ = i^„(0)/k„. As in the Szego theory [35] r„ and f„ are related to the above 
Toeplitz determinants by 



In[w{C)] 



rn = {-ly 



InHC)] ■ 



In |16| we have the correspondences of a{n) = (p„(0), a{n) = (pn{0), a{n) = k„-|_i/k„ and b{n) = rn- The 
Toeplitz determinants of central interest can then be recovered through the following result, which was given 
by Baxter in 1961 and is the generalisation of Eq. (11.16) of Geronimo and Case [TO] . 

Proposition 4.3 f[5].|15]). The sequence o/{/„[w]}J^q satisfy the recurrence 

In+l[w]In~l[w] 



(/,.H)^ 

subject to the condition r„f„ ^ 1 for all n > 1. 



= !-'"«?'„, n > 1, 



Further identities from the Szego theory that generalise are those that relate the leading coefficients back 
to the reflection coefficients. For example we have 



(4.9) 



<<5n(0)((5„(0), 



which is an extension of Eq. (11.10) in |16j . 

A fundamental consequence of the orthogonality conditions is a system of coupled linear first order dif- 
ference equations. They can constitute the starting point for the theory of a bi-orthogonal polynomial 
system rather than the orthogonality conditions and this was how Baxter [5] , [3] developed his theory. These 
recurrence relations are equivalent to Eq. (II. 7, 8) of Geronimo and Case [151 . 
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Proposition 4.4 ([21 The polynomial pair ^„(z),(^*(z) satisfy the coupled recurrence relations 

(4.10) K„V3„+l(z) = Kn+lZLpn{z) + (/3„+l (0)(^* (z) , 

(4.11) Kn^l^i{z) = Kn+lLpl{z) + l^„+i (0)z(^„ (z) . 

Together equations (|4.10I4.11I) arc equivalent to a single second order linear differenee equation and there- 
fore admit another linearly independent solution. We define another polynomial solution pair ■i/)„(z), (z) 
- the polynomial of the second kind or associated polynomial 

(4.12) -0„(z) := /" 7;^$^w(C)[Vii(C) - Vn(2)], n>l, f/jQ := kqWo = I/kq, 
and its reciprocal polynomial ^/;*(z). The integral formula for -0* is 

(4.13) C(^) := - / ;5^7^«^(C)[^"^«(C) - <(^)], " > 1, V^* := l/^o- 

A central object in our theory is the Carathcodory function - the generating function of the Toeplitz elements 

(4.14) F{z) := I ^}-±^w{0, ziT, 



It 27r< C- z 

which defines the inner and outer functions with expansions inside and outside the unit circle respectively 



F<{z)=wo + 2j2wkZ^ if|z|<l, 

fe=i 
00 

i^>(z) = -«;o-2^«;_fcZ-^ if |z| > 1. 



Having these definitions one requires two non-polynomial solutions eri(z), e* (z) to the recurrences (|4. 1014. lip 
and these are constructed as linear combinations of the polynomial solutions according to 

(4.15) e„(z) 0„(z) + Fiz)cp„{z) = / . ^ + ^ ^C)^„(C), 

(4.16) eUz) := ^:(z) - F(z)<^:(z) = - - / i^^wiO'fiUO, 

for n > 1 and eq = ^□(u'o + F), eg = ko{wo — F). Equation (|4.16l) along with the definition (|4.5p . is 
directly comparable to Eqs. (2.10), (2.14) and to Eqs. (2.11), (2.15) of the Wu study [32] in their respective 
z-regimcs. 

Proposition 4.5 f [17 ) .[T5].}19 | .[2S].|15|). The associated functions e„(z),e*(z) satisfy a variant of the cou- 
pled recurrence relations ^JTT^JTTl^ namely 

(4.17) K„e„+i(z) = K„+ize„(z) - .^„+i(0)e* (z), 

(4.18) K„e*+i(z) = K„+ie*(z) - (p„+i(0)ze„(z). 

The linear independence of the two solution sets to the coupled recurrences (|4. 1014. lip has the following 
consequences. 

Proposition 4.6 (|17|.[T5]). The Casoratians of the polynomial solutions ipn, </?^, ipn, V'n of the polynomial 
and non-polynomial solutions (p* , e„, e* are 

(4.19) ipn+l{z)lpniz) - 'lpn+liz)(Pn{z) = (/3„+l (z)e„ (z) - e„+i (z)(p„ (z) = 2 '^"+^ ^^-^ ^" ^ 

(4.20) ^:+,{zmz) - i'l+iizmz) = <+i(^)<(^) - <+i(^)<(^) = 2^^^^z"+\ 

(4.21) (^„(z)^;(z) -f V^„(z)(^;(z) = V?„(z)e:(z) + e„(z)(^;(z) = 2z". 

We note that the recurrence relations for the associated functions en{z), ^niz) given in (|4.17l4.18p dif- 
fer from those of the polynomial systems (|4. 1014. lip by a reversal of the signs of (pn{0), <^n(0). We can 
compensate for this by constructing the 2x1 vectors 

which can be directly compared with the polynomial solution 4'*(Z, 71) = Yn{z) in [16| . 
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Our extension of the matrix recurrence relations Eq. (11.12, 13) of Geronimo and Case [12] is the following 
result. 

Corollary 4.1 ([H]). The recurrence relations for a general system of bi-orthogonal polynomials Il4-10\4-n\ l 
and their associated functions ^.17\4-rS\ l are equivalent to the first order matrix recurrence 

(4.23) r„+i = A„y„ ^ f ^;+^(«))r„. 

A consequence of is that Kn has the property det JiTn ~ z. 

In Geronimo and Case |16| the matrix A*{n) — A'„. In addition to this matrix they defined the matrix 



*(Z, n) 



which figures in an alternative matrix recurrence of the matrix variable 

.^->;(^). 

Accordingly we can also define an alternative matrix of polynomial solutions as 

(4.24) *„(z) - ( "^"^^^ 
which satisfies a variant of (|4.23p . namely 

(4.25) ^,^^, = ^(_ -"+1^ ^n+i(0)z" 

The utility of this definition is that the determinant of the matrix given above is unity. 

The analogue of the Christoffel-Darboux summation formula is given by the following result. 

Proposition 4.7 ([3],[I1]). The summation identity 



K ■ ) 2^^]{z)ipj(i,) 
holds for 7^ 1 and n > 0. Here 

^(C) = C>n(l/C)- 

Another related bilinear identity is the formula for the discrete Wronskian. 

Proposition 4.8 ([16j). The discrete Wronskian of two solutions ^'i,n,^'2,)i to the recurrence system |^.g5[ ) 
is 

i.e. W[^i^n+i, ^2,n+i] = ^2,n] • Naturally this latter relation still holds for non-polynomial solutions 

of [rJE\ ) as well. 

Because we will be expanding the recurrence relation solutions about z = 0, oo we give the leading order 
terms in the expansions of (/9„(z), (^* (z), e„(z), e* (z) both inside and outside the unit circle. 

Corollary 4.2 f [15j). The bi-orthogonal polynomials (^„(z), (^* (z) have the expansions 

f(/9„(0) + ^(A^„^„_l(0)+^„(0)r„_l)z + O(z2) |z|<l, 

(4.27) ipn{z) = < Kn-1 

[k„z" + ;„z"-i + o(z"-2) kl > 1, 

K„+r„Z + 0(z2) |z| < 1, 

(4.28) ^:(z) = { _^^^p^^„ ^ _^(^^^^_^(o) + ^„(0)/„_i)z"-i + 0(z"-2) |z| > 1. 

Kn—1 

The associated functions have the expansions 

(4.29) !^,<(,) = ,«_i!^,"+i + 0(z"+2), |z|<l, 

Kn >, , (/5„+l(0) -1 , / </'«+2(0) (p„+i(0) /„+i 



(4.30) ^<{z) = ^'^^^^"' ^-^ + ^-^^y-' _ ^"^^v"^ --^^ ) z-' + 0{z-~'), |z|>l, 
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(4.31) ^<<(z) = + ( ^^I1±M _ ^j^h^) ,„+2 ^ o(^„+3)^ < 1^ 



(4.32) ^e:>(z) ^ 1 - hl+L;,-^ + ( _ , ^ ^ 



2 n \ / ' \ 

Kn+1 \Hn+2l^n+l 

In the notations of Wu [12] we make the identifications 

2U{z) = ~Sn,o - K„0-"e;<(z), 2y(z-i) = -2 + K„e;>(z). 

Two crucial relations are the jump conditions on the associated functions at \z\ = 1. 

Proposition 4.9. The associated functions defined in the interior of the unit circle ej^(z), e*^(z) differ from 
those in the exterior (z) , (z) on \z\ — 1 through the jump conditions 

(4.33) w{z)^^{z) = -ie>(z) + ie<(z), 

(4.34) w{z)^U^) ie:>(z) - ie:<(z). 

The latter relation (|4.34p is equivalent to Eq. (2.16) of [42] which states that 

c{Ox{C) = i + cu{c) + v{c'), 

on Id = 1. 

We make the observation on the weight (|1.6p with T > Tc or ^ < 1 that 

/I c-l -l\l/2 /I r 1 \ 1/2 /I ^ X 1/2 



Clearly w'^-^'^" (z) has non-zero winding number and so Szego's theorem is not applicable, however w has 
zero winding number and is given by the corresponding weight for T < Tc with ^ 1— > and z 1— )■ z"^. 
This latter transformation simply implies Wn ^ W-n and the Toeplitz determinant is unaffected. However 
the modification of the weight by a rational factor w{z) i— > z^^w(z) induces a Christoffel-Uvarov-Geronimus 
transformation. The elementary Christoffel-Uvarov-Geronimus transformation and it inverse w n- z^^w are 
known from [22], |40] - for w 1— >■ zw this is the K = l,L = 0,a— )-0 case in (40] 



/„ ^ /+ = (-l)"r„/, 



n 1 
_2 ^" 



'n+1 



2 

r,r rJ = r„ . 



-+ 1 

rn ^ - — 



whilst for w i-^. z this is the iiT^ 0, L=l,/3—!-0 case 

/n ^ = (-i)"r„/„, 



K„ 1-^ K„ , l,'*,! j — Kn~ : 

rn+1 

1 

1-^ r„ = — , 



__ _ _ l^n-l Tn+lTn-l 
fn ^ T,, — fn . 



5. Scattering States and Jost Functions 

We will make a number of assumptions, some of which we place on the symbol and some on reflection 
coefficients directly, in order to proceed with our analysis - 

(i) The symbol w{z) is absolutely integrable in the sense ■u;(e*^) G Lp[— 7r,7r] for p > 1, 

(ii) That the symbol has zero winding number and logw(z) satisfies the condition of Szego's theorem 

00 

^ bllcpp < 00, 
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(iii) Corresponding to Eq. (111.18) of [TB] we assume that 

oo 

\/\rnrn\ < oo. 

(iv) In addition we assume lim„_^oc k„ = Kqo so that 

lim = 1, 



and as a consequence we have 



One can check that these facts are verified in our application to the Ising diagonal correlations in an a 
posteriori demonstration, however we do not pursue those details here. 

Definition 5.1. Let h{z) have a Laurent expansion in z. We define the non-negative, positive and negative 
parts of h as 

oo oo — 1 

n—0 71—1 n— — oo 

Definition 5.2 ([IB]). Let us define the scattering states for n > by 

(5.1) ^-w-fl""!'!)' *-«(^~)-frfl 

\(l)+n{z)J \<P-n{z) 

They satisfy (|4.25p subject to the boundary conditions 

(5.2) lim |0+„-z"i =0, \z\ < 1, 

n— >-oo 

(5.3) lim |0+„| =0, \z\ < 1, 

n— ^-OD 

(5.4) lim |0_„| =0, |z| > 1, 

71— >00 

(5.5) lim !</>_„ - = 0, |z| > 1. 

n— f oo 

Their linear independence is guaranteed by \E'_|_„] = 1 for all n > 0. We have the exact correspon- 

dences with Geronimo and Case, 4>±.n{z) = 4>±{Z,n) and 4'±_n{z) — 4>±{Z,n). 

On |z| = 1 we decompose the polynomial solution 'i'n of (|4.25p into '^±n components 

introducing the Jost functions f±{z) as the coefficients of this decomposition. Inverting these relations one 
finds that 

f+{z) = -ipn{z)$+n{z) + (fn{z~^)(p+n{z), 
f-{z) = 'fin{z)(f>-n{z) - (fn{z^^)^-n{z)- 

As a consequence we have the limiting relations, which are a direct analog of Eq. (III. 13) in [IF] 

f+{z) = lim z"(^„(z-i) =: ^*^{z), \z\ < 1, 

n^oo 

/_(z) lim z-'V„(z), \z\ > 1. 

■n— >oo 

These limits can be shown to exist because /+ = VF[\['„, ^+„] and /_ = — M^[^'„, 4'-ri], and the property of 
the discrete Wronskian given in Proposition 14.81 We note the special values /+(0) = /_(oo) = Kqo- In the 
notations of Wu [32] we can identify P{z) = f+{z) and Q{z^^) = f-{z). In the notations of Borodin and 
Okounkov [6] we identify f+{z) = exp(— [y(— z)]>o) and f-{z) = exp(— [V^(— z)]<o). 
A central orthogonality result of [16] that generalises without change is the following. 

Corollary 5.1 ([IS]). Let us assume f+{z) ^ for \z\ < 1 and f-{z) ^ for \z\ > 1. For n > m the 

hi- orthogonal polynomials {(^„,(^„}5^q satisfy 

dz </5„(z)(^„i(z~i) 



T 27riz f+{z)f^{z) 
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= 0. 



Thus we have a factorisation formula, up to a constant, for the weight 
(5.6) w{z) cx 



f+iz)f-(.z) 

Proof. This follows from the proof given in [TB] with the obvious modifications. □ 

The scattering solutions defined in Definition 15.21 can in fact be related to the associated functions previ- 
ously defined in ^ 

Proposition 5.1. We can identify the scattering solutions i5.1\) together with the boundary conditions h5/A 
15.5)) with the associated functions 

(l)+n{z) = \f+{z)e<{z), \z\ < I, 
0+„(z) = -i/+(z)z-"<<(z), |z| <1, 
<^_„(z) = -i/_(z)e>(z), \z\>l, 
</>_„(z) = i/_(z)z-"e:>(z), \z\ > 1. 

Proof. Given that the combinations on the right-hand sides satisfy the relation (j4.25p the correspondences 
follow from a comparison of the boundary conditions (|4. 29114. 3^ with (j5. 2115. 51) . □ 

We will require expansions of the scattering states in their respective domains of analyticity, with respect 
to the monomials 



(l^+niz) = 


n' >n 




i+n{z) = 


^2(", 

n'>l 


n')z^', 


(j)-n{z) = 


n' >n 




4>-n{z) = 


Y ^2(": 


n')z-"' 



In conformance with the argument made in |16j we need to define a scattering function. 
Definition 5.3. The scattering function is defined on \z\ = 1 as 

and the sequence of Fourier coefficients Fm, Fm, m e Z by 

F^ := I ^z-Siz), F^ := / ^z" ' 



2'Kiz Jrp 2Tiiz S{z) 

In comparing with the work of Borodin and Okounkov [5] we have S{z) = exp(— y*(2)). For systems of 
orthogonal polynomials on the unit circle one has |>5'(z)| = 1 on |z| = 1 and Fm is the complex conjugate of 
Fm however this is no longer true for bi-orthogonal polynomial systems. 

We can now deduce a system of linear relations for the coefficients defined above which constitute a 
discrete analog of the Marchenko integral equations in the inverse scattering theory. 

Proposition 5.2. The coefficients Ai, A2, Ai, A2 satisfy the linear relations 



(5.7) 




= Ai{n, 


m) - 


h Y ^2('^, 

n'>l 




+rnj 


m > n > 1, 


(5.8) 





= A2{n, 


m) - 


h Y 


,n')Fn 


'+mi 


m,n > 0, 


(5.9) 




= Ai{n, 


m) - 


h Y ^2('^> 
n'>l 


n')F^. 


+rni 


m > n > 1, 


(5.10) 





= A2{n, 


m) - 


h Y 


.n')F„ 


'+mi 


m, n > 0. 



n'>n 
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Proof. These relations follow from a recasting of the jump equations (|4.33l4.34p using the factorisation result 
of CoroUarv lS.ll and Proposition l5.ll Thus to deduce (|5.7II5.10|) we employ 



(5.11) 
(5.12) 
(5.13) 
(5.14) 

respectively. 



— — - = (j)+n[z) + ——(/)_„ Z , 

f-iz) f-{z) 

^ = 0_„(.) + -^0+„(z), 



-^=0-„(.) + ^0+„(.), 



-^=0+„(z) + -^0_„(z) 



□ 



We further define kernel matrices Gi^m,Gi^m (which differ from those given by Eq. (V.ll) of [IS] being 
the transpose of theirs - this appears to be a typographical error) using the scattering Fourier coefficients 

1.' tr 



7n' > 1 

We are now in a position to solve the linear system given in the previous proposition. 



m'>l 
7n' > 1 



Proposition 5.3 (|16]). The solutions of the discrete Marchenko equations ( [5. 7^5. lO\) for the norms and 
reflection coefficients are given by 



(5.15) 

where 



det[l + G]^ = det[l + G]: 



In addition we have 



dot 



_ det[l + G]- 1 
nl det[l + G]- ' 



/1 + G 

Gn+l,n 1 + Gn+l,n+l 



G 



n+2,n 



G. 



Gri.ri+2 Gn.n+d, 
Gn+lji+2 G„+i^„+3 
1 + G„+2,Ti+2 G„+2,ri+3 



■A 



(5.16) 



det[l + G]- 1 



(5.17) 



'n+l 



det[l + G]-i 





^ Fn+l 


Fn+2 


Fn+3 


...\ 




Gn+l,n 


1 + Gn+l,n+l 


Gn+lji+2 




det 


Gn+2,n 


Gn+2,n+l 


1 + G„+2,ji+2 






Gn+3,n 


Gn+3,n+l 


Gn+3,n+2 












■•/ 




^ Fn+l 


Fn+2 


Fn+3 


...\ 




Gn+l,n 


1 + Gn+l,n+l 


Gn+l,n+2 




det 


Gn+2,n 


Gn+2,n+l 


1 + G„+2,n+2 






Gn+3,n 


G„+3^„+l 


Gri+3.n+2 






\ 









Proof. For convenience we define the coefficient ratios 

Ai(n, m) 



a{n, m) 



Ai{n,n) 



a{n, m) 



Ai{n, m) 
Aiin, n) 
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Summarising (j5.7ll5.1(J)) we now have for ?7i > ?i > the Unear relations 

oo 

(5.18) ^ a{n,m) + G,n.n + ^ Gm.ja{n,l), 

oo 

(5.19) ^ a{n,m) + Gm.n + ^ Gm.iHnJ), 

l=n+l 

and for m = n > 

oo 

(5.20) + E G^MnJ), 

oo 

(5.21) ^ ^^1 + G„.„+ V Gn.ia{n,l)- 

From the solutions for Ai, A2, Ai, A2 one can recover the leading polynomial normalisation and reflection 
coefRcients 

—2^ = Ai{n, n) = Ai{n, n), 

A2{n,l) _ A2{n,l) 
Aiin^n) Ai[n,n) 

Using Cramer's rules we can solve the linear relations (|5.18l5.19p for a(n, Z), a{n, I) in terms of determinants 
and the results are (|5. 1515. 1615. 17p . □ 

The Neumann expansion of the basic determinant 

det[l + G]^ ~ 1 + 51 Gn,,n, + det[G„^,„,],,fc=i,2 + . . . , 

ni>n n2>ni>n 

is directly applicable to the formulae given in SjSl (|3.4p and p.5p . after noting that Gij = —K{i,j). It is this 
object which specifies the Toeplitz determinant (|4.ip 

_ lo det[l + G]- 



det[l + G]if 



6. Equality of the Fredholm Determinants 

6.1. A-Extended Riemann-Hilbert problem. We now pose a new problem which is an extension of the 
jump conditions given in Proposition 14.91 bv introducing a parameter A through the mappings ej^ 1— > A^e^ 
and e*^ i— > A^e*^. We have two new functional equations where the associated functions defined in the 
interior of the unit circle e^(z), 6*^(2:) differ from those in the exterior e^(z),e*^(z) on \z\ = 1 through the 
jump conditions 

(6.1) w{z)ip^{z) = -ie>(z) + A2ie<(z), 

(6.2) w{z)^l{z)^\e*>{z)^\^lel<{z). 

Given a weight 'w{z) such a system of functional equations can be taken as the defining relations of the 
bi-orthogonal polynomial system and this is essentially the Riemann-Hilbert approach (see Proposition 2.8 
of [TS]). The solutions of this extended Riemann-Hilbert problem, , (/s* , e„ , e* , acquire a A dependence as 
do their coefRcients and therefore the reflection coefficients. However we do not address here the question 
of a A dependent weight from which one could construct these coefRcients - this will be the subject of our 
companion work |41| . 

Having set ourselves the task of solving this new problem we retrace the steps taken in Section [5] that 
logically follow from using Equations (|6.ip and (|6.2p . We find that Proposition l5.2l now becomes the following. 
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Proposition 6.1. The X-extended coejjicients Ai, A2, Ai, A2 satisfy the linear relations 

— 5ra,n = Ai{n,m) + ^ A2{n,n')Fn'+mi m > 71 > 1, 



= A2(n, m) + ^ ri')K'+m, m, n > 0, 

n' >n 

— ^m,n = m) + ^ A2(n,n')-Fn'+m, m > n > 1, 

n' >1 

= ^2(71,771) + ^ Ai(n, n')Fn/+,„, 771, 71 > 0. 



Proof. The alteration to the proof of Proposition [521 consists of tlie replacement of (|5.12p and (|5.13p by the 
equations 

-^=0-„(.) + A-^0,„(z), 
respectively and the same reasoning applied. □ 



Consequently Proposition 15 . 31 is extended in the following way. 

Proposition 6.2. The solutions of the X-extended discrete Marchenko equations for the norms and reflection 
coefficients are given by \5.15\) . \5.16\) . {5.11^ respectively with the substitutions Fm,F„i 1— > XFm,XF„i o,nd 
Gm,mGm,n ^ A^Gm,„, A^Gm,„. As a conscqucnce the coefficients have expansions 

^2 

(6.3) ^ 1 — A Gn.n + A > Gn.niGni,n + . . . , 

" Tli>Tl 



(6.4) rn+i ^ Ai^„+i — A'^ > F„j+iG„j „ + . . . , 

A— !-0 ^ — ' 

ni >n+l 

(6.5) r-n+l ^ AFfi+i — A'^ > Fni+lGni,n + . ■ . ■ 

A-i-O ^ — ' 

ni >n+l 

Proof. For 777, > 77 > the linear relations (|5.18p and (|5.19p are replaced by 

00 

= 0(77,771) + A^G„i,„ + A^ ^ G„i,ia{n,l), 

l=n+l 

00 

= 0(77,777) + A^Gm,n + A^ ^ Gm,ia{nJ), 

l=n+l 

respectively and for 777 = 77 > Equations (|5.20p and (|5.2ip are replaced by 

00 

j4 r = 1 + A^G„,„ + A^ Gn,ia{n,l), 

A^rvn -ill [Ti, Tl) 

00 

= 1 + A2G„,„ + A2 GnMnJ), 



K^Ai{n,n) ^^^^^ 
respectively. These are solved in the same way as indicated in the proof of Proposition 15. 31 □ 
The Neumann expansion of the basic determinant is now 

det[l + A^G],^ 1 + E + E det[G„^,„J,- fc=i,2 + • . ■ , 

ni>n n2>ni>n 



and this specifies the extended Toeplitz determinant (j4.ip 

^ Iq det[l + A'G]^ 
" ~ det[l + A2G]S°' 
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6.2. Iterative solution of A-extended jump conditions. In the final step of the proof of Proposition [3TT] 
we return to the expression of the form factors as multiple integrals given in fJTl The calculation of the form 
factors p.l4ll.l6p by Lyberg and McCoy follows the method of Wu |12 which is essentially the application 
of Wiener-Hopf arguments combined with an iterative solution to the jump condition (|6.2[) . The last part 
of our proof formalises this assertion. 

Corollary 6.1. Consider the low temperature case ^ > 1. With 

Uiz) = (1 - r'z)'/', f-{z) = (1 - kr'z-'r^/^ 

the solution of the coupled set of equations 

(6.6) K^^Uz)^f+{z) ([/_(z)]>„-A2iK„ [/-W<<W]>o) , 

(6.7) i«nC(^) - 1 = (- [f-{z)U,+X''^Kr. [/-(^)<<(^)]<o) > 

(6.8) «:„z>:(z-i) =/_(z-i) ([z"/+(^-')]>o+ [^"/+(^-')(^^n<>(^^')- l)]>o) : 

(6.9) i«„z"e:<(z-i) = j-^^ ([^"/+(^-')]<„ + [z"/+(z-i)(i^„4>(z-i) - 1)] J , 

developed as an expansion in gives the form factor expansion formula il.l3\) with 1^1. 

Proof. Employing the usual Wiener-Hopf arguments we can deduce the coupled functional equations 

(6.10) K^^Uz) = f+{z) {[f-{z)U - ¥n [/-(^)<<(^)]>o) , 

(6.11) \^ne:>{z) 1 = (- [f-{z)U + l^n [f-{zV<{z)\^^) , 

(6.12) «.„z>:(z-i) =/_(z-i) ([z"/+(^-')]>o+ [^"/+(^-')(^^«4>(^-')- l)]>o) : 

(6.13) i«„z"e:<(z-i) = j-^^ ([z"/+(^-')]<„ + [z-f+{z~'){\n^e*>{z-') - 1)] J , 

from and These are precisely Equations (2.19a), (2.20a), (2.19b) and (2.20b) of [H] respectively, 

when re- written in our notation. The first, second and fourth equations can be considered as a coupled set of 
functional equations for (^* (z), e*^(z) and e*^(z). We now seek an iterative solution to the second and fourth 
equations starting with the initial term involving e^{z) neglected in (j6.1ip . i.e. a perturbation expansion 
with the parameter as a book-keeping parameter recording the order of approximation in this solution. 
This means that we replace (|6.10p and (|6.1ip with (|6.6p and (|6.7p respectively, which is expressed by the 
mapping e*^(z) i— ^ A^e*^(z). This is the essential idea behind the asymptotic analysis employed by Wu 
and also the derivation of the form factor expansion by Lyberg and McCoy, who employed Cauchy integral 
representations for the positive and negative parts, and which ultimately led to the formulae (|1.14p and 
(|1.16p . By introducing A in this way we can trace it through the workings of |28| and find that it appears 
precisely as the correct pre-factor of the terms in the form factor expansion. □ 

This completes the proof of Proposition 13. II and we close this part of our study by making two remarks. 

Remark 6.1. The two Fredholm determinant evaluations for the diagonal correlations of the Ising model 
reported here don't have direct relevance to the multi-variable integral formulas for the form factors of 
correlations away from the diagonal, such as those in [43) or |32j . and more work needs to be done to relate 
these apparently different representations. 

Remark 6.2. Both our Fredholm determinant formulae generalise beyond the diagonal correlations of the 
Ising model and are valid for symbols or weights that are of regular semi-classical type. Therefore all of our 
results would apply to the row correlations as well with the obvious modifications. 

6.3. Explicit examples. The general theory given above can now be applied to give explicit formulae for 
the Ising model diagonal correlations in the low temperature regime ^ > 1, due to the fact that the winding 
number vanishes and the conditions on the weight and scattering data ensure the validity of this theory. 
We will evaluate a number of initial data explicitly and this allows us to make contact with the results in a 
companion study [41] where the same quantities have been calculated using the fact that our bi-orthogonal 
polynomial system is also an isomonodromic system identified with Picard's solution of the sixth Painleve 
equation. 
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According to Corollary 15.11 for | ^| > 1 we factorise the weight (|1.6p into the interior and exterior Jost 
functions 



U{z) = (1 - Uz) = (1 - i-'z-'y 



1/2 



In this example Kqo = 1 as cq = 0. A straight-forward calculation reveals that the Fourier coefficients with 
1^1 > 1 and n G Z are 

^" = ^^'rl^t^ rl"l2fi(i \n\ + i; \n\ + 1; r^), 
TT ' |n|! 

ZTT ' |n|! 

All these coefficients have complete elliptic integral evaluations (using the standard definitions given in 
Chapter 19 of [31] or Chapter 13 of [12]) as the following examples illustrate 



(6.14) Fo = -K{t), 

TT 

(6.15) Fi = F_i = [K{t) - E(0] , 

TT 

(6.16) Fo = -[(t-l)K(i) + 2E(0], 

TT 

(6.17) A = F_i = -^t-'^' [{t - l)K(t) + (< + l)E(i) 
The preceding evaluations imply that the kernel matrix l,m G Z is given by 

:-2„r(n + /-i)r(n + m+i) 



l,m. 



1 ,.-2» r(n + /- 2 



27r ^ (n + l)l(n + m)\ 

71 — 1 ^ ^ ^ 

X 2-Fl(-i n + Z - i;n + ? + l;<)2Fl(i 77, + 777 + i; 77 + 777 + 1; t), 



which is convergent for |^| > 1 as 

F„ ^ ^-1/2(1 -i)-l/2|„|-l/VH/2, F„ ^ _ 1^-1/2(1 „t)l/2|„|-3/2^-|"l/2, 
n— J-±oo n— fiboo 

However because of the general summation identity Z, 777 G Z, ? 7^ 777. |t| < 1 

E°° r(77 + ^ — ■^)r(77 + 777 + ■^) , , , , 

t" , , I' , ^2i^i -i?^ + ;-5;r7 + / + l;i 2i^l in + 777 +i; 77 + 777 + l;i 
_ (77 + i)!(77 + 777)! z / z ^ 



~2tt 



("l)i+l(^)m+l ^ 

?!777! Z — 777 



^ 2i^i(i^ + i;^ + l;t)2i^i(im+ i; 777 + 2; i) 



777 + 1 

1 



l + l 



2J 



Fi (i , 777 + I ; 777 + 1; t)2Fi (i , Z + i ; Z + 2; t) 



we find Gi^m ~ ^K{1, m), as given by (j3.3p . This provides an alternative derivation of the discrete Fredholm 
determinant, to that given in [5]. The diagonal elements can be computed from the recurrence relation 



(6.19) Gij - i-n[±iHi±ilt'+i2i^^(_ 1 , ^ + 1 ; Z + 2; t)2F,{^,l + |; ? + 2; i) + Gz+i,,+i, 



and the identity Go.o + = —1. 
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Again we find that the kernel matrix has complete elliptic integral evaluations as the early examples 



1 n n ipfifp 


\f\ < 1 






[b.ZU) 




[tt + 8b(r)K(i) + 4(i — l)K(t) J , 




(6.21) 




-1 = — r ' — 41?j(r)K(r) — z(i - 


- l)K(t) J , 


(6.22) 


n 1 r 


-TT^ + 4(t - l)K(t)^ + 8K(t)E(t)J , 




(6.23) 


G*— 1,1 = — SGi. 


_i = ^<-i [E(i) - K(t)] [(t + l)E(i) + 


(t-l)K(t)], 


(6.24) 




= r 6E(t)2 _ 4(t _ 2)E(t)K(t) 


+ 2(i-l)K(t)2] , 


(6.25) 




^ [-STT^t + 4(t - l)(3t - 2)K{tf + 8(3t 


- 2)K(t)E(t) + 8(< + l)E(i)2' 



Remark 6.3. An important feature of the kernel matrix, given by (|6.18p and (|6.19p . that emerges is the 
appearance of the regularised Gauss hypcrgeometric function rather than the function itself and the conse- 
quence is that the matrix exists for negative indices /, m, i.e. for all /, m € Z. This in turn means that it is 
possible to give definite meaning to the normalisation coefficients k„ and the reflection coefficients r„, f„ for 
negative indices. This is an important theme in our companion study |41| . 

These examples furnish a check on our results and a link to exact evaluations of the bi-orthogonal polyno- 
mial system coefficients using the connection with Picard's solution to the sixth Painlevc equation, see [29] . 
|41j . From these works it is known that 

(6.26) ^ = — = sec(x) [cn(z, t)dn(z, t) + sn(z, t)Z{z, t)] , 

in the low temperature phase where the independent variables are 

2K(0 . . , , 
z := X, A = sm(a;), 

and A = 1 corresponds to x = 7r/2. The cn(z, t), dn(z, t), sn(z, t) are the standard Jacobian elliptic functions 
and Jacobi's elliptic zeta function is 

n 

with Jacobi's elliptic epsilon function defined by 

E{z,t) := E(am(z,t),<), 

where E{z,t) is the incomplete second elliptic integral and the elliptic amplitude function am(z,t) is defined 
through the inversion of the incomplete second elliptic integral 

z := F(am(z, t), t). 

By expanding the exact result (|6.26p about A = wc find that 

t " ^ + i " ~ " 8E(t)K(t)] A^ 

-t- 2^ [Qtt'' - 407r2(t - l)K{tf + l&it^ + 2t- ■i)K{tf + E(i)(-807r2K(t) + 64(t + l)K{tf)\ A^ 

+ 0(A«). 

Clearly this coincides with (|6.3p for n = and (|6.22p to O(A^). In addition it is known from the same works 
that 

(6.27) ^ = = sec(a;) ^ 



cn(z, t)dn(z, t) + sn(z, t)Z(z, t) ' 
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and the leading terms in the expansion of this ratio give that 



77 ^ ^ + i + 8E(i)K(0 + 4(t - l)K(t)2] 

+ 0(A'^), 



247r4 



which coincides with (|6.3p for n = — 1 and (|6.20p . 

Other important examples are ro,fo in the low temperature phase, which arc no longer unity for A 7^ 1 
but rather [IT] 



(6.28) 
(6.29) 



1 



ro = 



ro = sn(z,i), 

[1 + cn{z, t)dn(z, t) + sn{z, t)Z{z, t)] [1 — cn(z, t)dn{z, t) — sn(z, t)Z{z, t)] . 



sn(z, t) 

In the limit A — > 1^ we note that ro,fo 1 using the elliptic functions expansions 



cn{z,t) = x/l~lK(f)(a;- f) + 0(x- f)3, 



dn(z, t) = + 0{x - f )^ 

sn(z,t) = l + 0{x- f)2, 



2 

TT 



Z{x, t) = -- [E{t) + {t-l) K{t)] (x - f ) + 0{x - f )^ 



By expanding (|6.28p and (|6.29p about A = we have 



-K(i)A + -i^K(t) [tt' - Ait + iMtf] A3 + OiX% 



and 



ro 



1 



;2(t-l)K(t) + 4E(t)]A 



+ — [-24E(t)2K(t) + 2E(t) (tt^ - 12(t - l)K{t)^) + {t - l)K(i) (tt^ - 4(t - l)K(t)2)] A^ + 0{X^). 

These coincide with (|6.4p and (|6.5p respectively to 0(A) using the evaluations (|6.14p and (|6.16p . 

We can also make contact with the form factor expansions with the Appell function kernel given in fj^l 
and without any loss of generality we examine the case n = 0. From the works of Orrick et al [35], Boukraa 
et al [7] and Mangazeev and Guttmann [29] we have expressions for the non-trivial boundary values Iq for 
A^l 



(6.30) 



[l-t) 



1/4 



MO 


T)'dl{x 


r) 


MO 


r)?94(0 


r) 



t^e, ^ < 1 



where we adopt the standard definitions of the elliptic theta functions (see Chapter 20 of [3T] and §13.19 of 
[12j ) and the elliptic nome is defined by 



and T = 



Naturally in the case A = 1 we have /q = 1 as one can verify from (|6.30p using the identities 



t = 



^2 



MO\t) 



l-t = 



MO\t) 



MO\t) 



, MhMr)=MO\r), MhT^\r)=MO\r). 



Expanding the expressions (|6.30p about A = we have in the low temperature regime 
(1 - t)-'/"/o = 1 + A^{t) [K{t) - E{t)] X' 



^K(i) - ^E(t) + ^Kitf + iK(i)E(t)2 - ltK(t)3 - iK(t)2E(<) 



X^ + OiX"^), 
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which concurs with the initial terms of (|2.8p with (|2.10p and (|2.9|) . In the high temperature regime we have 



{l-t)-'/^Io = -K{t)\x+ — 



^-i(i-2)K(t)2-iK(t)E(i) 



16 



3^^ 



12 



— - -K(t)E(t) - -(i - 2)K(t)2 ( - - mm ) + ^mfm' + ^ {t 



120 



6t + 6) K{tY 



which agrees with the initial terms given by (j2.13p with (|2.14[ 12.151 I2.16p . The form factor expansion 
coefficients of this correlation function have been given in §4 and §5.2 of [31], and we observe agreement with 
the leading terms given above. In addition these coefficients have been given in Appendix A of |29| . 
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